Quantum phase transitions in an interacting atom-molecule boson model 
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We study the quantum phase transitions of a model that describes the interconversion of inter- 
acting bosonic atoms and molecules. Using a classical analysis, we identify a threshold coupling line 
separating a molecular phase and a mixed phase. Through studies of the energy gap, von Neumann 
entanglement entropy, and fidelity, we give evidence that this line is associated to a boundary line 
in the ground-state phase diagram of the quantum system. 
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I. INTRODUCTION 



A Quantum Phase Transition (QPT) is characterized 
by a change in the properties of the ground state of a 
system (i.e. at zero temperature) as some parameter is 
varied across a critical point [l|,y]- This parameter can be 
an external magnetic field as in the quantum Hall effect 
or superconducting materials or, as in a Bose-Einstein 
condensate, a change in the the s-wave scattering am- 
plitude or external fields to produce an atom-molecule 
condensate 3] . The change in the nature of the ground 
state is typically identified by non-analyticity of some 
quantity, such as the ground-state energy or a correla- 
tion function. A common approach for studying QPTs is 
to use concepts such as order parameters and symmetry 
breaking a la Landau-Ginzburg, in analogy with thermal 
phase transitions. However in many cases of interest or- 
der parameters are difficult to identify. More recently, 
in the wake of quantum information theory, it has been 
realized that one can employ alternative concepts such as 
entanglement y and fidelity [3,0 as means of identifying 
different phases. 

Bosonic models are useful testbeds for investigating 
QPTs in that they accommodate large particle numbers 
with few degrees of freedom, in contrast to fermionic 
systems, which simplifies the analysis. In this paper 
we study QPTs for an interacting atom-molecule bo- 
son model. The model does not possess any symmetries 
which might give rise to symmetry-breaking order. Us- 
ing other methods, we will establish a line of QPTs in 
the ground-state phase diagram. This line is first iden- 
tified via phase space bifurcations in a classical analysis 
of the systenL, adapting a correspondence that has been 
put forth in Q relating bifurcations in classical systems 
with entanglement entropy of quantum systems. We then 
confirm that this line is associated with quantum phase 
transitions of the quantum system through studies of the 
energy gap, entanglement entropy, and fidelity. Even 
though a QPT is only rigorously defined in the thermody- 
namic limit iV — > cx), where N is the number or particles 
in the system, we will show that the aforementioned con- 
cepts do respond to changes in the ground state of the 



system for finite TV and strongly indicate the presence of 
QPTs (cf. i). 

We consider the model for interacting atomic and 
molecular bosons described by the following Hamiltonian 



i 



H 



UaNl + UbNl 



UabNaNb + fiaNa + fibNb 



-fl{a^a^b + b^aa). 



(1) 



where a^ is the creation operator for an atomic mode and 
b^ is the creation operator for a molecular mode. The 
parameters fii are chemical potentials for species i and 
fl is the amplitude for the interconversion of atoms and 
molecules. The Hamiltonian commutes with the total 
atom number N — Na + 2N},, where Na = a^a and Nj, = 
b^b. Thus the Hamiltonian can be block diagonalized into 
sectors on which N takes a constant value. Hereafter we 
consider the action of the Hamiltonian to be restricted to 
such a sector and treat A^ to be a scalar variable. Notice 
that the change of variable il ^ — f7 is equivalent to the 
unitary transformation 

b -^ -b. (2) 

The parameters Uj describe s-wave scattering, taking 
into account the atom-atom (Ua), atom-molecule (JJat) 
and molecule- molecule {Ub) interactions. 

In the no-scattering limit Ua — Uab — Ub = 0, the 
model (P) has been studied using a variety of methods 
[l3-[l3|i and it was found to undergo a quantum phase 
transition when Ha/^ = v2N. This was confirmed ana- 
lytically through a Bethe Ansatz solution [lj|. However 
in the experimental context the s-wave scattering inter- 
actions play a significant role |15l4l9j . It was shown in 
|20l - |22| that for the general model ([T]) the inclusion of 
these scattering terms has non-trivial consequences for 
the physical behaviour. In this paper we re-examine the 
QPT found in the no-scattering limit, showing how the 
concepts of quantum information, entanglement and fi- 
delity are related to it. Moreover, we also investigate the 
effect of the s-wave scattering parameters in the QPTs, 
establishing a connection between the QPTs and the bi- 
furcation line associated to the global minimum of the 



Hamiltonian ([T]) in the parameter space of the classical 
system. 

The paper is organized as follows: in section 2 we give 
an outline of the classical analysis of the model ([1]) includ- 
ing the fixed point bifurcations that occur in the classical 
phase space. In section 3 we investigate the behaviour 
of the energy gap, entanglement, and fidelity to iden- 
tify QPTs. We also establish a connection between these 
QPTs and a bifurcation line of the global minimum in the 
parameter space of the classical system. Our conclusions 
are stated in section 4. 



2- 




II. CLASSICAL ANALYSIS 



We first motivate undertaking a classical analysis to 
gain insights into the existence of QPTs in quantum sys- 
tems. Recall the proposition in Q that, if for a bipartite 
quantum system there exists a supercritical pitchfork bi- 
furcation of the global minimum at some coupling param- 
eter in the phase space of the analogous classical system, 
then the quantum system will take a maximum value 
of the ground-state entanglement entropy at that cou- 
pling. This proposition holds true for attractive bosons 
in a double- well potential, within a two-mode approxima- 
tion. The double-well model admits a QPT associated 
with symmetry-breaking as shown in |23| . The critical 
coupling coincides with the supercritical bifurcation of 
the global minimum in phase space [2J], and also coin- 
cides with the point where the entanglement entropy is 
maximal [25|. 

To investigate the extent to which these ideas extend to 
the present model we briefly recall the classical analysis of 



the Hamiltonian ([T|) as given in !2Q.| . Let Nj 



j ^a,b 



be quantum variables satisfying canonical commutation 
relations. We make a change of variables from the oper- 
ators j, j^, j ~ a, b via 



j = exp(i6'j)^iV,-, f = ^Nj cxpi-iOj) 



such that the canonical commutation relations are pre- 
served. Now define the variables 



z^-iN,-2Nb), 



N 
= -i29a-9t), 



FIG. 1. Parameter space diagram identifying the different 
types of solution for equation ((9|. In region I there are no 
solutions for z when 6 = 0, and one solution for z when 9 = 
Nn/4:. In region II there are two solutions for z when ^ = 0, 
and one solution for z when 9 = Nn/A. In region III there 
exists one solution for z when ^ = 0, one solution for z when 
9 = Nir/i, and two solutions for 9 when z = —1. In region 
IV there is one solution for z when 9 = 0, and no solution for 
z when 6 = Nn/A. In this region the global minima occur at 
z = —1 and all values of 9. The boundary separating regions 
I and II of the region III is given by A = a -|- 1, while the 
equation A = a — 1 separates the region III of the region IV. 
The boundary between regions I and II has been obtained 
numerically. 



where 




(4) 
(5) 
(6) 



We note that the unitary transformation ^ is equiva- 
lent to6 —>■ B+Nir/A. Hereafter we restrict our attentions 
to the "repulsive" case A > 0. 

We now regard ([3]) as a classical Hamiltonian and in- 
vestigate the global minima of the system. The first step 
is to find Hamilton's equations of motion which yields 



dz dH 4 ^-— 



de 

'di 



dH 
dz 



2Xz + 2a 



-z)il 
l-3z 



^)sin(f ), (7) 



(8) 



such that (z, 9) are canonically conjugate variables. In 
the classical limit where N is large, but still finite, we 
may equivalently consider the transformed Hamiltonian 

m 



A9 



H = Xz^ + 2az + /3 + ^2(1 - z)(l + z) cos ( — ) (3) 



Within the interior of the compact phase space the fixed 
points of the system are determined by the condition 



dz d6 

-T = -T = 0. 
dt dt 



(9) 



Extremal points may also occur on the boundaries 
z = I, z = — 1. The bifurcations of the fixed points di- 
vide the coupling parameter space into different regions. 






FIG. 2. (Color online) Dimensionless energy gap between 
the first excited state and the ground state as a function of 
a = fia/^V'^N for different values of A^ and A = 0. Here 
Q.^1, fib = and Ua = Ub/-i = 0.25. 



as shown in Fig. [TJ Of these bifurcations, only the bound- 
ary separating regions III and IV (the hue A = a — 1) is 
associated with a quahtative change of the global mini- 
mum of the Hamiltonian ([3]) . In the regions I, II, III qual- 
itative changes are associated with either saddle points 
or maxima (see [20| for details). For these regions the 
minimum of ([3]) occurs at = Ntt/A with z > — 1. In 
region IV there is a line of global minima for (jH]) at the 
boundary z = —1 and for all values of 0. These global 
minima do not satisfy ^. 

From the above we see that there is a bifurcation of 
the phase space minimum from region III, where there 
is a unique minimum, to region IV, where there is a line 
of minima. The bifurcation line in the parameter space, 
which is given by A = a — 1, is not of a supercritical 
pitchfork type discussed in [3]. Consequently, we should 
not expect the ground-state entanglement entropy to be 
maximal on this bifurcation line. In fact it has been 
shown in [TlJ that for the point (A, a) = (0, 1) on this line 
the entanglement entropy is not maximal. Nonetheless, 
we will establish below that the bifurcation line A = a — 
1 is still associated to a line of QPTs for the quantum 
Hamiltonian ([TJ. 



III. QUANTUM PHASE TRANSITIONS 



A. Energy gap 



We begin our analysis by considering the dimension- 
less energy gap between the first excited state and the 
ground state, ^Ej [ig.- Using numerical diagonalization 
of the Hamiltonian, in Fig. [2] we plot the scaled gap 
^E I [iai as a function of the coupling a, for A = and 
different values of N . We observe that as iV increases the 
dimensionless gap decreases and the coupling approaches 



FIG. 3. (Color online) Dimensionless energy gap between the 
first excited state and the ground state as a function of a = 
^,a/^\/2N for different values of A and iV = 2000. Values of 
Q,, Hb and Ua, Ub are the same as in the previous figure. These 
results indicate that the minimal values lie approximately on 
the line A = a — 1. 
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FIG. 4. Ground-state expectation values for the atomic num- 
ber fraction Na/N and the molecular number fraction Nb/N 
as a function of q = ^a/^V2N for A = and TV = 2000. We 
have set f2 = 1, ^b = and Ua ~ Ub/4: — 0.25. There is a 
sharp transition at a = 1. 



the value a = 1. Fig. |3] shows similar results for fixed 
N = 2000 and varying A. We observe that the occurrence 
of the minima of the gap, determining the QPT, fit well 
with the predicted boundary separating regions III and 
IV given by A = a — 1 . This is the first evidence suggest- 
ing that a line of QPTs occurs. Note that in both graphs 
the scaled gap goes to zero at a single point, rather than 
over an interval, of the coupling a. This is indicative of 
the fact that there is no phase where the ground state is 
degenerate, which would be a requirement for the exis- 
tence of a broken-symmetry phase. 

To better understand the physical meaning of the 
QPTs, we depict in Fig. |4] the ground-state expecta- 
tion value of the scaled atomic number operator (solid 
line) and the scaled molecular number operator (dashed 
line) as function of a for iV = 2000 and A = 0. The aver- 
age value for the number of atoms decreases while the 
average number of molecules increases as a increases. 
For a > 1, the average number of molecules is maxi- 
mal, therefore we can interpret this point a = 1 as the 
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FIG. 5. (Color online) Entropy of entanglement of the ground 
state as a function of a = /la/^V^N for different values of N 
and A = 0. Here n = 1, ^6 = and Ua = t/b/4 = 0.25. 
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FIG. 6. (Color online) Entropy of entanglement of the ground 
state as a function of a = /la/^V^N for different values of A 
and N = 2000. Here n = 1, pi, = and Ua = Ub/A = 0.25. 



threshold coupling for the formation of a predominantly 
molecular state. This result is consistent with the clas- 
sical analysis, whereby in region IV the global minima 
have z = —\ which corresponds to a molecular phase. 



B. Entanglement 

One may consider the atom-molecule model ([1]) as a 
bipartite system of two modes, A and B. In this case, the 
standard measure of entanglement is the von Neumann 
entropy of the reduced density operator of either of the 
modes [Il| . The state of each mode is characterized by its 
occupation number. Using the fact that the total number 
of atoms N is constant, a general state of the system can 
be written for even iV in terms of the Fock states by 



N 



\^>)=Y,dn\2n)\N-n) 

n=0 

where d„ are complex numbers. 



(10) 



FIG. 7. (Color online) First (insert) and second derivative of 
the entropy of entanglement of the ground state as a function 
of Q = iia/9.y/2N for iV = 1000 and A = 0. Here fi = 1, 
/ii, = and Ua = Ub/4, = 0.25. 



The density operator for state pU)) is given by 

N 

p=\^){^= Y. d*^d^\2m)\N-m){2n\{N-n\. (11) 

m,n— 

Taking the partial trace with respect to the mode B 
yields the reduced density operator for the mode A, 



N 



PA 



Tr, 



Kp) = E 



|2n)(2n| 



(12) 



n=0 



The entropy of entanglement of the ground-state of the 
system is given by 



N 

E{pa) = -Tr[pA\og^{pA)] ^~Y. \dn?^og^(\dn?). 

n=0 

(13) 
Using the above expression (fT5|) and the energy levels ob- 
tained through exact diagonalization of the Hamiltonian 
^ we plot in Figs. ([5]) and (jH]) the entropy of entangle- 
ment of the ground-state as a function of the coupling a 
for different values of a and N. 

In Fig. ([5|) we confirm that the entanglement entropy 
is not maximal at the threshold coupling (A, a) = (0, 1) 
determined from the classical analysis, in agreement with 
[ll| . However we do observe that the entanglement en- 
tropy exhibits a sudden decrease close to a = 1 that be- 
comes more pronounced as the total number of atoms 
N increases. Fig. ([6]) shows similar results for fixed 
N — 2000 with varying A. In this latter case we see 
that the occurrence of the abrupt decay of the entropy 
fits with the predicted boundary separating regions III 
and IV given by A == a — 1. 

We also depict in Fig. ([7]) the first and second deriva- 
tives of the ground state entanglement entropy as a func- 
tion of a. Combined, Figs. ([5]) and ([7]) suggest a discon- 
tinuous behaviour in the limit iV — > oo, consistent with 
the existence of a QPT. 




FIG. 8. (Color online) Fidelity of the ground state as a func- 
tion of Q = ^ia/^V^N for n = 1, A = 0, iV = 1000 and 
Ua ~ Ub/4: = 0.25. We keep ^t ~ for the first state and 
vary ^t — ^ for the second state. 
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FIG. 9. (Color online) Fidelity of the ground state as a func- 
tion of Q = fia/^V^N for A = 0, ^i, = (first state) and 
fib = 1 (second state) with varying A'^. Here = 1 and 
Ua = Ub/4: — 0.25. The insert shows the minima close the 
critical point a = 1. 



C. Fidelity 



Another possibility to investigate QPTs is through the 
behavior of the fidelity [3, Q. This concept is widely 
used in quantum information theory [20|. The fidelity 
is defined as the modulus of the wavefunction overlap 
between two states 

j-(Va0) = |(VI0)I- 

In Fig. ([5]) we present the wavefunction overlap between 
two ground-states corresponding to the external param- 
eter Hb = for one ground state and fib = 1 for the other 
ground state. For 7 = the states are indistinguish- 
able and there is no information about the QPT. When 
7 increases the states become distinguishable and the fi- 
delity has a minimum at the same point where we find 
a QPT using the energy gap and entanglement entropy. 
For larger values of 7 the distinguishability increases but 
the minimum value remains in the same position. 

Fig. dni) depicts the ground state fidelity \{fib = 0|^6 = 
1)1 for A = 0, 7 = 1 and varying N. With increasing 
N the states become more distinguishable and the point 
where the minimum occurs moves towards a = \. Fig. 
pn)) shows similar results for fixed N = 2000 and varying 
A. In this case we also observe that the occurrence of the 
minima of the fidelity, determining the QPT, fit well with 
the predicted boundary line A = a — 1. 



IV. SUMMARY 
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FIG. 10. (Color online) Fidelity of the ground state as a 
function of q = fj,a/^V2N for iV = 2000, fib = (first state) 
and fib ~ i (second state) and varying A. Here Q — 1 and 
Ua = Ub/4: = 0.25. 



tations in the limit N — > cx3. Also, the derivatives of 
the entanglement entropy rapidly vary in the vicinity of 
this line. Finally, the fidelity approach shows that the 
states with fia = and fia = j become distinguishable 
on this line, where fidelity has a minimum. The param- 
eter 7 changes the distinguishability of the states but 
does not change the value of the threshold coupling. All 
approaches indicate markedly different behaviours when 
the line A = a — 1 is crossed, giving a strong indicator 
towards the existence of a line of QPTs. These results 
show that the study of bifurcations in the phase space of 
classical systems to identify quantum phases transitions 
may be applicable at a general level, independent of the 
nature of the bifurcations. 



We have used three different approaches to show that 
the interacting atom-molecule boson model ([1]) exhibits 
a line of QPTs that coincides with the boundary of the 
parameter space diagram determined by bifurcation line 
A = Of — 1 of the global minimum of the Hamiltonian 
in the classical analysis. First, the dimensionless energy 
gap is minimal on this line and indicates gapless exci- 
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